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PROCEDE D« ESTIMATION DU MOUVEMENT 
ENTRE DEUX IMAGES 

5 La presente invention concerne un precede d ' estimation 
du mouvement entre deux images numeriques. 

Le mouvement entre deux images success ives, Ii et I2, 
se definit generalement sous la forme d»un champ de 

10 mouvement associe I'une ou 1' autre des images I^, I2 
et constitue de vecteurs deplacement se rapportant 
chacun a un point de 1 ' image concernee- Le vecteur 
deplacement est un vecteur a deux dimensions 
representatif de la difference de position entre le 

15 pixel de 1 ■ image Ii et le pixel associe de 1» image I2 
se rapportant au meitle point physique de la scene 
filmSe. 

Une estimation de mouvement est utile dans les domaines 
20 du traitement de 1' image n^cessitant la connaissance 
des mouvements ou des disparitSs entre deux images. On 
peut citer, a titre d'exemples, les domaines 
d* applications suivants: 

- le codage video : le champ de mouvement defini & 
25 partir d' images deja decodees est alors utilise pour la 

prediction de 1 • image suivante; 

- 1« imager ie medicale : le procede d* estimation du 
mouvement entre deux images est utilise pour faire une 
analyse de mouvement du coeur par exemple; 

30 - la telesurveillance : le procede peut servir au 
controle du trafic routier; 



- la reconstruction tridimensionnelle a partir d' images 
laulti-vues : le precede est utilise pour estimer les 
disparitSs entre les differentes vues. 

Pour obtenir un tel champ de mouvement, il est connu de 
decomposer 1 • image en elements finis. Ces elements 
finis qui peuvent etre par exemple des triangles ou des 
quadrangles, sont determines par un maillage dont les 
noeuds correspondent aux sommets des elements finis. Un 
vecteur deplacement est calcule pour chaque noeud du 
maillage. Puis, par le biais d'une fonction 
d» interpolation, on peut en deduire un vecteur 
deplacement pour chaque point de 1 » image consideree. Le 
champ de mouvement est ainsi determine par un modele 
d'elements finis, definissant le maillage utilise pour 
partitionner 1* image en elements finis et la fonction 
d» interpolation permettant de calculer le vecteur 
deplacement en tout point de 1' image. 

Le maillage utilise peut etre regulier ou non et doit 
etre choisi suffisamment dense pour model iser au mieux 
le mouvement entre les deux images sans toutef ois 
necessiter une quantity de calculs ou de donnees a 
transmettre trop importante. Ce choix est effectue une 
fois pour toutes en debut de precede et ce maillage 
reste generalement le meme tout au long de 
1» estimation. 

Le calcul des vecteurs deplacement des noeuds du 
maillage peut §tre effectue selon diverses methodes. II 
existe tout d'abord les methodes dites de mise en 
correspondance, ou "matching", consistant a tester un 
ensemble discret de valeurs possibles de vecteurs 
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deplacement pour chaque noeud du maillage et a retenir 
les meilleurs vecteurs selon un critere donne. Une 
seconde methode dite par transformee consiste a 
utiliser les proprietes de la transformee de Fourier et 
5 de ses extensions af in de convert ir le mouvement en un 
dephasage dans l*espace transforme. Enfin, il existe 
une troisieme methode dite dif f erentielle permettant de 
determiner les vecteurs deplacement par optimisation 
d'un critere mathematique (par exemple, une erreur 

10 quadratique entre 1 • image et sa valeur predite avec l-e 
champ de mouvement) . Cette methode est la plus utilisee 
pour 1* estimation de mouvement avec mod^lisation par 
elements finis. Une methode dif f erentielle classique 
pour 1 'optimisation des vecteurs deplacement est la 

15 methode de Gauss-Newton. La presente demande concerne 
plus particulierement la famille de precede 
d* estimation de mouvement utilisant un modSle 
d* elements finis et une methode dif f erentielle pour 
deiterminer le champ de mouvement. 

20 

Bien que tres repandu, ce type de precede presente 
cependant quelques inconvenients. Le maillage choisi en 
debut de procede peut se reveler inappropriS par 
rapport au contenu semantique de 1' image, certaines 

25 zones de 1' image meritant un maillage plus dense et 
d'autres un maillage plus aere. Par ailleurs, sous 
I'effet du champ des vecteurs deplacement des noeuds du 
maillage, le maillage initial sur 1' image de depart, 
par exemple 12, se transforme en un nouveau maillage 

30 sur I'autre image, par exemple II. 11 peut alors se 
produire des situations patholdgiques au niveau du 
nouveau maillage telles que : 



- des retournements d' elements finis: des elements 
finis se retournent et viennent en recouvrir d'autres, 
detruisant ainsi la propriete de partition du domaine 
de 1' image que doit verifier tout maillage, 

- debordement des noeuds peripheriques du maillage 
deplacfe apres application des vecteurs deplacement au- 
deia du domaine de 1 ' image II: certains pixels de 
!• image 12 peuvent etre associes a des pixels de 
!• image II situes en dehors du domaine de 1' image 11^. 
Ceci n'est pas forcement genant mais il peut etre 
interessant de forcer les noeuds peripheriques du 
maillage deplace a rester dans le domaine de 1' image 
II. 

Un but de 1 ' invention est de proposer un procedS 
d' estimation de mouvement dans lequel le maillage est 
optimise au cours de 1' estimation de maniere a obtenir 
en fin de proced§ un maillage adapte au contenu 
semantique des images. A cet effet, on prSvoit un 
affinage local des elements finis durant 1" estimation 
de mouvement, 

Un autre but de 1' invention est d'ameliorer 
l»efficacit§ de la methode de Gauss Newton pour 
1« optimisation des vecteurs deplacement des noeuds du 
maillage, A cet effet, on prevoit d'effectuer cette 
optimisation sur plusieurs niveaux de resolution des 
images • 

Enfin, un autre but de 1' invention est de proposer un 
precede d' estimation de mouvement permettant d'eviter 
les situations pathologiques precitees* A cet effet, on 
prevoit selon 1» invention de raj outer, au cours de 



l"etape d ' optimisation des vecteurs deplacement, des 
contraintes pour eviter de telles situations. 

Aussi, 1" invention a pour objet un precede d' estimation 
du mouvement entre deux images numeriques, et I2, de 
luminance et Y2, destine a generer, pour chaque 

point de coordonnees x,y de 1 ' image 1 2 un vecteur 
deplacement d(x,y)=(dx,dy) de maniere a former une image I2 
a partir de !• image Ii, de luminance Y2(x,y)=Yi(x-dx,y-dy}, 
qui soit une approximation de 1' image I2/ 
caracterisg en ce qu'il comporte les Stapes suivantes: 

(a) - definir un modele d» elements finis initial 
comportant un maillage dont les noeuds sont des points 
de 1' image I2, un vecteur deplacement a chaque noeud 
dudit maillage, et une formule d' interpolation pour 
calculer la valeur du vecteur deplacement de chaque 
point de 1" image I2 a partir des valeurs des vecteurs 
deplacements des noeuds de la maille a laquelle il 
appartient, 

(b) - optimiser la valeur des vecteurs deplacement du 
modele selon une methode dif f 6rentielle, 

(c) - calculer un ecart E entre 1 ' image ±2 et 1 • image 
I2 pour chaque element fini ou maille, 

(d) - operer un maillage plus fin sur une fraction 
discrete de 1» ensemble des elements finis 'determinee 
selon un critere relatif aux ecarts E, et attribuer un 
vecteur deplacement a chaque nouveau noeud die maillage, 

(e) - rep§ter les etapes (b) , (c) et (d) sur le modele 
d' Elements finis obtenu au terme de I'etape (d) 
prScedente, jusqu'a satisfaire un critere d'arret. 
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Selon un mode dc realisation ameliore, on definit en 
outre, pour chaque image nuitierique Ii et I2/ un 
ensemble de R images If de niveau de resolution r et de 

luminance Y{ avec r prenant les valeurs (0, /R-1) et 

i les valeurs 1 et 2, les images I? et I? correspondant 
aux images num§riques Ii et I2, et on effectue les 
etapes (b) a (e) pour chaque niveau de resolution r, a 
partir du niveau r=R-l jusqu'au niveau r=0. 



Enfin, selon un mode de realisation prefere, on prevoit 
raj outer des contraintes sur le deplacement des 
elements finis au moment de 1 • optimisation des vecteurs 
deplacement de maniere a §viter le retournement des 
15 Pigments finis, Selon un autre mode de realisation, on 
prSvoit egalement d'introduire des contraintes pour 
6viter le debordement du maillage obtenu apres 
application des vecteurs deplacement au-dela du domaine 
de 1» image II, 



D'autres caracteristiques et avantages de 1 ' invention 
apparaltront a la lecture de la description d6taill§e 
qui suit et qui est faite en reference aux dessins 
annexes, dans lesquels: 

- la figure 1 represente un schema d'un premier 
mode de realisation du precede d' estimation de 
mouvement selon 1 " invention; 

- la figure 2 illustre I'etape (d) du procedS de 

1 • invent ion ; et 

- la figure 3 represente un schema d'un mode de 
realisation ameliore du precede d" estimation de 
mouvement selon 1' invention. 



On considere deux images numeriques Ii et I2 de 
luiainance respective Yi et Y2. proc6d6 de 

1» invention consiste a generer, pour chaque point P de 
coordonnees (x,y) dans !• image I2, un vecteur 
deplacement d(x,y)={dx,dy). Ce vecteur est defini conime 
fetant le vecteur permettant de construire a partir de 
!• image Ii, une image ±2 de luminance Y2(x,y)=Yi(x-dx,y-dy) 
qui soit une approximation de I2. Les deplacements sont 
ainsi dSfinis de 1 • image Ii vers 1* image I2. 

Le champ de mouvement recherche est defini par un 
module d"elements finis. Dans la suite de I'expbse, on 
considere que les 616ments finis sont des triangles 
sans qu'on puisse y voir une limitation de la portee de 
la prfisente demande S cette forme d' elements finis. Le 
modele d' elements finis comporte par consequent un 
maillage triangulaire, des vecteurs deplacement definis 
aux noeuds du maillage, lesquels noeuds correspondant 
aux sommets des triangles, et une formule 
d' interpolation pour calculer le vecteur deplacement 
des points a I'intferieur des triangles. 

La formule d' interpolation utilisfee pour calculer le 
champ de mouvement en tout point du domaine de 1 ' image 
I2 est la suivante: 

si on considere le point P de coordonnees (x,y) dans 
1» image I2 appartenant au triangle e de sommets P^, Pj 
et Pk de coordonnees respectives (Xi,yi), (Xj,yj) et (xk,yk), son 
vecteur deplacement est egal a 

d(x,y)=: i;xi/f(x,y)-d(x,,yi) 

oil reprSsente une fonction de base associfee au 

triangle e. 
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Dans le cas d'une interpolation affine, les vf(x,y) 
repr^sentent les coordonnSes barycentriques du point P 
dans le triangle e de sommets Pi, Pj, Pjc- Ces. fonctions 
sont d^finies par la relation suivante : 

v|/i(x,y)=ai+Pix + yiy (x.y)ee 

l=i.j,k 



,M/f (x.y) = 0 



(x.y) €e 



Xjyfc - x„yj -t- (y i - yic )x + (x^ - Xj)y 
^'^"''^=xjy.-x.yi^x,y,-x,y.H-x,y,-x,y, 

Les fonctions affines Vj(x,y) et Vk(x,y) se diduisent de 
la fonction M/f(x,y) en permutant circulairement les 
indices i,j,k. On peut ^galement utiliser des modules 
d'614ments finis plus €volu§s, les fonctions v pouvant 
alors fetre €tendues a des polynomes de degr^ n>2 mais 
la formule d' interpolation des vecteurs dSplacement 
fait alors intervenir des dSrivSes premiferes, 
secondes. .. Un choix divers de modeles d 'Pigments finis 
est proposS dans I'ouvrage ""-Handbook of Numerical 
Analysis" de P.G Ciarlet et J.L Lions, Volume 2, pp 
59-99, Ed. North Holland. 

Selon 1' invention, au fur et a mesure de 1 • estimation 
de mouvement, on optimise la valeur des vecteurs 
dSplacement des noeuds du maillage appeles vecteurs 
nodaux, et on densif ie localement le maillage lorsque 
c'est n^cessaire. Avantageusement, cette optimisation 
sera effectuee sur plusieurs niveaux de resolution en 
partant d'un niveau de resolution faible. 



Selon un premier mode de realisation illustre par la 
figure 1, le procedfe de !• invention comporte cinq 
etapes r6f erencees (a) a (e) . 

Selon I'etape (a), on definit un modele d'61ements 
finis initial en selectionnant des points de 1" image I2 
selon un maillage triangulaire initial . Les noeuds du 
maillage representent les sommets des triangles 
(Elements finis) du modele, Ce maillage peut Stre 
d§fini de maniere quelconque, par exemple en fonction 
des besoins de 1 ' application ou de connaissances 
prealables ou du mouvement deja calcule entre deux 
images prScedentes de la meme sequence video. Si aucune 
donnee concernant le maillage initial n'est specif iee, 
on utilise un maillage regulier en quinconce. Les 
mailles sont alors des triangles. Un vecteur 
dSplacement de valeur nulle est alors associe a chaque 
noeud du maillage. La formule d • interpolation dSfinie 
prScedemment est Sgalement une donnfee du modfele 
initial. 

Selon I'etape (b) , on optimise la valeur des vecteurs 
deplacement du modele selon une mSthode dif f erentielle, 
par exemple la methode de Gauss-Newton ou son extension 
de Marquardt. Cette optimisation peut etre soit libre^ 
c*est-a-dire sans contraintes imposees sur les valeurs 
possibles des vecteurs nodaux, soit avec contraintes. 
Les vecteurs nodaux designent les vecteurs deplacement 
des noeuds du maillage. L" optimisation avec contraintes 
est directement liee a 1 ' optimisation libre et fait 
I'objet d'un mode de r§alisation developpe plus loin. 
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La technique d» optimisation libra employee ici exploits 
les caracteristiques interessantes de la m^thode de 
Gauss-Newton (convergence rapide de 1 •optimum) et de la 
mfethode du gradient avec pas adaptatif (convergence 
globale vers un optimum local) pour rSsoudre le systeme 
lin6aire a suivre. Cette technique est une correction 
iterative des vecteurs deplacement d(x,y) permettant 
d'obtenir des le depart une approximation grossiere du 
mouvement. Le nombre d* iterations k P^nr cette 
optimisation des vecteurs deplacement est soit spScifie 
par I'utilisateur en debut de procSd^, soit fpnction 
d'un seuil li€ a la variation maximale entre deux 
valeurs consecutives de vecteurs nodaux pour deux 
iterations successives. Nous developpons ci-apres 
1» extension Marquardt de la methode d ■ optimisation de 
Gauss Newton. 

L'expression des corrections SD^""* a 1' iteration k+l 
des vecteurs deplacement selon cette methode est donnSe 
par le systeme systeme lineaire suivant: 

avec : 

- D*"*^ une matrice colonne de 2N elements comprenant les 
composantes dx et dy des vecteurs nodaux a 1» iteration 
k+l, N §tant le nombre de noeuds du maillage a I'etape 
courante ; 

- D*" une matrice colonne de 2N elements comprenant les 
composantes dx et dy des vecteurs nodaux a 1» iteration 
k; 

- H = [R^+a.l2N] 



la matrice identite de dimension 2N; 



matrice colonne de 2N elements dont 

N elements V^^nE^ et N Elements Vy , n 
dSsignant un noeud du maillage et 
prenant tour a tour les valeurs (1...N); 
avec 



\ „ „ 9Ii(x-dv,y-dv) ^ 

Vx.„E''=2. X E DFD,(x.y). ^ .M/'n(x.y) 

cesupp(n)(x,y)ec 



ax 



C6supp(n)(x,y)ee 



5y 



oCl DFDk(x.y) = Y2(x,y)-Y,(x-d>,,y-dy) a 1 ' iteration k 
oCl supp(n) represente le support de la f onction de 
base Vn(^»y) attachee au noeud n, c'est-a-dire 
1» ensemble des triangles ayant le noeud n pour 
sommet; 



une matrice carrSe de dimension 2N 



Rmr=2- 5: Z ( 

e € siipp(inn)(x.y) ee 



eesupp(mn) (x,y)€e 

ail(x-dj^,y-dy) ai|(x-d^,y-dy) 



ax 



ay 



).H/^{x.y).M/n(x.y) 



pk,yx _ T> k,xy 
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eesupp(inn) (x,y)ec 

oa m et n dSsignent des noeuds du maillage et 
prennent tour a tour les valeurs (1...N) et oil 
supp(nin)=supp(n)o supp(m). 

- a = inax„(|V„E''|.||xi/nl|) 

oH Iv/nl est une norme f onctionnelle de \|/n • Les 
deux normes les plus int^ressantes sont: 



IVnIl = SUP(x.y)«up p(n)kn(x,y)| = 1 OU 

Ikni = j| Z[Vn(x,y)f 

" \|supp(n)|(x.y)esupp(n) 



10 |supp(n)| dSsigne le cardinal de la region discrete 

supp(n). 

A la fin de cette phase d • optimisation, on dispose de N 
vecteurs nodaux se rapportant chacun a un noeud du 
maillage. 

15 

Selon un mode de realisation en variante du gradient 
adaptatif, on peut 4galement envisager d'utiliser une 
technique de decomposition dite "LDL*" profil" en 
langage technique pour rSsoudre le syst4me lin^aire (1) 
20 et acceierer le traitement. Cette technique est dfecrite 
dans I'ouvrage intitule "Analyse num§rique matricielle 
appliquee a I'art de I'ingenieur" de Theodore Lascaux, 
Tome 1, pp 295-299, Ed- Masson,1986. 

25 Selon une caract€ristique importante de 1* invention, le 
maillage est ensuite affine localement par division des 
triangles lorsque I'Scart entre 1' image ±2 et 1' image 
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l2 sur ces triangles est trop eleve. Cest pourquoi, 
salon I'etape (c) du precede, on calcule un 6cart E 
entre 1' image I2 et 1' image I2 pour chaque triangle e. 
L'^cart E est defini de la maniere suivante 

E = 2]DFD-(x,y) 
(x.y)€e 

avec DFD(x,y) = Y2(x,y)-Y|(x-d^,y-dy) 

Bien entendu, pour calculer cet ecart pour chaque 
triangle, il faut avoir prealablement calcule la valeur 
des vecteurs deplacement de tous les points de 1' image 
I2 par interpolation a partir des vecteurs nodaux 
obtenus au terme de I'etape (b) . 

Ensuite, conformement a I'etape (d) , on affine le 
maillage sur une fraction discrete de 1' ensemble des 
triangles du modele, Cette fraction est dfeterminSe 
selon un critere relatif aux ecarts E calculus 
precedemment a I'etape (c) . Pour operer cet affinage, 
on peut par exemple classer les triangles du modele par 
ordre decroissant de leurs ecarts E et subdiviser les X 
premiers triangles de ce classement en triangles de 
taille inferieure. X est une fraction predeterminSe du 
nombre d' elements finis dans le modele, par exemple la 
moiti§. 

Pour af finer localement le maillage, on peut egalement 
prevoir de comparer 1' ensemble des ecarts E calcules a 
I'gtape (c) a un 6cart seuil qui est fonction de la 
taille de 1' element fini considere, et de subdiviser en 
elements finis de taille inferieure les elements finis 
dont les hearts E sont super ieurs a 1« ecart seuil. 
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La subdivision d'un triangle e en quatre triangles de 
taille inftrieure est illustre a la figure 2. Le 
triangle e est d^fini par les trois soinmets Pi, P2 et 
P3 ayant pour vecteurs deplacement respectif s dj , 62, 
dj. Pour le subdiviser en quatre, on d^finit trois 
nouveaux noeuds P4, P5 et Pe au milieu des trois c5t#s 
P1P3/ triangle. On attribue a chacun de 

ces trois nouveaux noeuds un vecteur deplacement §gal S 
la moyenne des vecteurs deplacement des deux soinmets du 
c6t€ auquel il appartient, respect ivement (dj+d3)/2, 
(dj + d2)/2, (d2+d3)/2. Les triangles adjacents au 
triangle e ayant pour c5te P1P2, P2P3 P1P3 sont alors 
eux-memes subdivisSs en deux ou trois. 

on obtient ainsi un module d«616ments finis dont le 
maillage a etS affinS localeihent. Selon I'etape (e) , on 
repete ensuite les Stapes (b) , (c) et (d) en prenant 
comme modele S optimiser le modele obtenu au terme de 
I'gtape (d) precedente. Et, on rSpete cette succession 
d' Stapes jusqu'a satisfaire un critere d' arret. Ce 
critere d'arrSt est par exemple un nombre predetermine 
d« elements finis a atteindre au terme de l»Stape (d) . 

on peut egalement prevoir d'arreter le procSde lorsque 
les ecarts E de 1* ensemble des elements finis du module 
obtenu au terme de I'etape (c) prScedente sont 
infSrieurs ^ un ecart seuil. 

Selon un mode de realisation amSliorS illustre figure 
3, on effectue les Stapes (b) a (e) en s • appuyant sur 
plusieurs niveaux de resolution d' images Ii et I2. 
cette variante a pour but d'amSliorer et d'accSlSrer la 
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convergence des calculs des vecteurs deplacement. Pour 
ce faire, on definit tout d'abord pour chaque couple 
d' images numeriques Ii et I2, un ensemble de R images 
I[ de niveau de resolution r et de luminance Yj^ , r 
prenant tour a tour les valeurs (R-l,R-2 , . • • 0) et i les 
valeurs (1,2) puis on effectue les etapes (b) a (e) 
pour chaque niveau de resolution r, a partir du niveau 
de resolution r=R-l jusqu'au niveau r=0. A noter que 
les images I? et I? correspondent aux images numeriques 
Ii et I2. 

En pratique, les images I- sont obtenues par filtrage 
de 1' image Ij au moyen d'un filtre passe-bas lineaire 
ne laissant passer que 1/2^ de la bande spectrale de 
!• image consideree dans les directions x et y, c"est-a- 
dire un filtre de reponse impulsionnelle ayant une 

bande passante BP3:=[-l/2^+l , 1/2^^+^] dans I'espace des 
frequences normalisees [-1/2,1/2]. L'image I[ est 
definie par la relation suivante: 

MM 

Yi'(x.y)= X 2:Yi(x-u.y-v)hX 

u=-Mv=-M 

Le filtre utilise est par exemple une approximation 

d'un filtre ideal et sa reponse impulsionnelle est 

definie de la maniere suivante: 

r M 
hl=^ avec S= -M<n<M 

S n=-M 



et 



; = 2B.sinc(27tBrn) = 2B— — ^ 



2' 

oa B et M sont des entiers naturels 
M =+00 dans le cas ideal. 
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Comme indique precedeiament , cette optimisation sur 
plusieurs niveaux de resolution permet d'ameliorer et 
d'acc^lSrer la convergence des calculs des vecteurs 
dSplacement. A noter que le nombre de niveaux de 
resolution R choisi peut diff^rer du nombre d'affinages 
successifs effectu^s sur le maillage. 

Selon un mode de realisation pr^ferS, on propose de 
rajouter des contraintes de compacite sur chaque 
triangle du modele de maniere a empScher le 
retournement des triangles. 

La compacit§ d'un triangle de sommets Pi, Pj, Pk est 
d6finie par la relation suivantei 

4*7t*S(Pi.Pj.Pk) 

C(Pi,Pj.Pk) - p2(p. p. p^) 

avec C(Pi,Pj,Pk)e ]0,1[; et 

S(Pi , P j , Pk) et P(Pi , P j , Pk) representant 
respectivement la surface et la perimStre du 
triangle (Pi,Pj,Pk)i 

Si on empeche la compacite d'un triangle de tendre vers 
z^ro, on I'empSche ainsi de se retourner. C'est 
pourquoi, pour 4viter les retournements de triangles, 
chaque triangle doit verifier la contrainte suivante : 

C(Pi + dp^ ,Pj + dp. ,Pk + dp^ ) ^ K *C(Pi.Pj.Pu) 

o K * C(Pi.Pj.Pk) - C(Pi + dp^ .Pj + dp. ,Pk + dp^ ) < 0 

o ge(D) < 0;e = triangle(Pi,Pj,Pk) 
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oa K est un param&tre f ixant la variation de compacitS 
autoris6e et D est le vecteur colonne des vecteiirs 
d^placenent des noeuds du modele. 

Selon 1' invention, on associe une telle contrainte a 
chaque triangle au moment de 1 • optimisation des 
vecteurs deplacement. L'etape d • optimisation des 
vecteurs deplacement revient a un systeme du type: 

ming E(D) 
^ ge(D)^0 VcgI 

oil -E(D) represente I'ecart entre 1' image I2 et 
1» image prSdite 

- ge est une contrainte liee au triangle e; 

- I est 1' ensemble des triangles du maillage. 

Pour resoudre les problfemes d' optimisation sous 
contraihtes, on privoit d'utiliser la technique dite du 
Lagrangien augmente. Cette technique est decrite dans 
I'ouvrage intitule '^^Theories et algorithmes » » de 
Michel Minoux, Tome 1, pp 257-260, Dunod 1983. Cette 
technique combine deux techniques d • optimisation: 
1" optimisation Lagrangienne et 1 • optimisation de 
p^nalites exterieures. 

Selon cette technique, resoudre le systfeme precedent 
revient a resoudre le systeme sans contraintes suivant: 

minf\ [E(D) + XG(gc(D)Ac.re)] 

oH re est un element de penalite, 

est un multiplicateur de Lagrange 



18 



G est un Lagrangien augmente determinS par la 
relation: 

\gc(P)+rSc(P) site >Oetge(D)SO 

Xege(D) sire=Oetge(D)>0 

0 si ge(D) ^ 0 



G(ge(D)Ac.re) = 



Les contraintes g© ont ete pr^alablement lin6aris€es 
par la formule de Taylor a I'ordre 1 : 

gc(D)«gc(Pi.Pj.Pk)+ Z dp. - (P.) 

Pi=Pi.Pj.Pk ^(Pi.Pj.Pk) 

La methode d' optimisation est alors la suivante:- 

- on initialise k=0 

- on pose ;i =0 et r=0, X e Ot"" et r e ^i™ , et m d^signe le 
nombre de triangles du modele 

- on determine ensuite le minimum 5D''**(A.,r) tel que, a 
1' iteration k+1, 

oil y^=(X ,r) , C* est une matrice de 91^^X91^"' 

C^Y forme une matrice des contraintes linear is6es 
ayant pour coefficients les valeurs suivantes: 



0 si G^f D) < 0 

aj g.<P,.Pj.P,, siP,=Pi.PjouP, 

Q sinon 



oH e reprSsente le triangle de sommets Pi, P j , Pjc. 

Puis, on effectue la mise a jour de X par I'algorithme 
d'Uzawa et on augmente r puis on repete 1' operation 
pr6c§dente jusqu'a ce que toutes les contraintes soient 
vSrifiees. On passe ensuite a 1' iteration k+2. A noter 
que I'algorithme d'Uzawa est developpe en detail dans 
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I'ouvrage intitule ""Theories et algorithmes • ' de 
Michel Minoux, Tome 1, Dunod 1983. 

Selon un dernier mode de realisation, on peut ^galement 
5 pr^voir d'introduire des contraintes pour eviter le 
dgbordement du maillage obtenu apres application des 
vecteurs deplacement au-dela du domaine de !• image II. 
Ce mode consiste a forcer les noeuds periphSriques du 
maillage a rester sur les bords de 1 " image aprds 

10 application des vecteurs deplacement. Pour ceia, les 
composantes abscisses 5D^^^ pour les noeuds 
p6riph€riques sur les bords gauche et droit del' image 
12 sont annulees a chaque iteration k. De meme, les 
composantes ordonnees SD*^"^^ pour les bords haut et bas 

15 de !• image 12 sont annulees a chaque iteration k. 



REVENDICATIONS 



1. Proc§d6 d' estimation du mouvement entre deux 
images numferiques, Ii et I2, luminance Yi et ¥3, 

destinS a g§n6rer, pour chaque point de coordonnSes x,y 
de 1' image I2 un vecteur deplacement d(x,y)=(d.x.dy) de 
mani&re a former une image a partir de 1 • image Ii, 
de luminance Y2(x.y)=Y,(x-d.x,y-dy). qui soit une 
approximation de 1' image l2r 

caract^rise en ce qu'il comporte les etapes suivantes: 

(a) - d^finir un modele d' elements finis initial 
comportant un maillage dont les noeuds sont des points 
de 1" image I2, un vecteur deplacement a chaque noeud 
dudit maillage, et une formule d- interpolation pour 
calculer la valeur du vecteur deplacement de chaque 
point de 1' image I2 a partir des valeurs des vecteurs 
dfeplacements des noeuds de la maille a laquelle il 
appartient, 

(b) - optimiser la valeur des vecteurs deplacement du 
modele selon une methode dif f erentielle , 

(c) - calculer un 6cart E entre 1 • image la et 1' image 
I2 pour chaque element fini ou maille, 

(d) - opferer un maillage plus fin sur une fraction 
discrete de 1' ensemble des elements finis dStermin^e 
selon un critere relatif aux Scarts E, et attribuer un 
vecteur dSplacement a chaque nouveau noeud de maillage, 

(e) - repeter les fetapes (b) , (c) et (d) sur le mod&le 
d' elements finis obtenu au terme de I'etape (d) 
precedente, jusqu'a satisfaire un critere d' arret. 
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2. ProcedS selon la revendication 1, caract§ris6 en ce 
que, pour operer un maillage plus fin sur une fraction 
discrete de 1» ensemble des 616inents finis a l»6tape 

5 (d) , on classe ledit ensemble des elements finis par 
ordre decroissant de leurs ecarts E et on subdivise les 
X premiers elements finis de ce classeroent en elements 
finis de taille inferieure^ X representant une fraction 
predfeterminee du nombre des elements finis de 
10 1» ensemble. 

3. Proc^dS selon la revendication 1, caract§ris§ en ce 
que, pour op§rer un maillage plus fin sur une fraction 
discrete de 1» ensemble des elements finis a I'etape 

15 (d) / on compare 1' ensemble des ecarts E calcules a 
I'etape (c) a un ecart seuil qui est fonction de la 
taille de 1' Pigment fini consider^, et on subdivise en 
felfements finis de taille infSrieure les Elements finis 
dont les hearts E sont super ieurs a l»6cart seuil. 

20 

4. ProcedS selon I'une des revendications 1 a 3, 
caracteris€ en ce que ledit critere d» arret est un 
nombre predetermine d' Elements finis constituant le 
modele d' Elements finis, a atteindre au terme de 

25 I'etape (d) . 

5. ProcSde selon I'une des revendications 1 h 3, 
caract6ris6 en ce que ledit critere d' arret de 1 ' fetape 
(e) est satisfait lorsque les ecarts E de 1' ensemble 

30 des elements finis du modele obtenu au terme de I'etape 
(d) precedente sont inferieurs a un ecart seuil 
fonction de la taille de 1' elements fini considere. 



6 ProcSde selon I'une des revendications 1^5, 
clract6ris6 en ce qu'on d^finit en outre, pour chaque 
image numerique H et I2, un ensemble de R images I[ de 
niveau de resolution r et de luminance Y;' avec r 

prenant les valeurs (0 R-D et i les valeurs 1 et 

2, les images I? et I? correspondant aux images 
numSriques Ii et l2» 

et en ce qu'on effectue les etapes (b) a (e) pour 
Chaque niveau de resolution r. a partir du niveau r=R-l 
jusqu'au niveau r=0. 

7 ProcedS selon la revendication 6, caract^risS en ce 
que les ensembles de R images de niveau de resolution r 
sont obtenus par filtrage des images II et 12 suivant 
les deux directions x et y au moyen d'un filtre passe- 
bas de rfeponse impulsionnelle hj, , chaque image If €tant 
dSfinie par la relation suivante: 

X'(x.y)= Z i:Yi(x-u.y-v)hX 

u=-Mv=-M 

avec M entier naturel. 

8 Proc6d6 selon la revendication 7, caracteris6 en ce 
que la reponse impulsionnelle K dfefinie de la 

maniSre suivante: 

r _,sin2KBrn 
h[, = ^ avec si, = 2B.sinc(2nB,n) = 2B-^^^ 

M 

S= Z^n 

n=-M 

B un entier naturel. 



9. Precede selon I'une des revendications precedentes, 
caracterise en ce que les vecteurs deplacement initiaux 
sont des vecteurs nuls. 

10. Proc6d6 selon l*une des revendications precedentes, 
caracterisfe en ce que I'ecart E entre 1' image ±2 et 
1' image I2 pour chaque 616ment fini e est defini par la 
relation suivante: 

E - XDFD-(x,y) 
(x.y)€c 

oQ DFD(x,y) = Y2(x,y)-Yi(x-d^,y-dy) 

11. Proc^de selon I'une des revendications prec§dentes, 
caracterise en ce que la formule d • interpolation pour 
calculer la valeur du vecteur deplacement d'un point P 
de coordonnSes (x,y) dans 1 ' image I2 appartenant a 
1« element fini e de sommets Pi, Pj et P^ de coordonnees 
respectives (x^yO, (Xj,yj) et (xk,yk) est la suivante: 

d(x,y)= XM'i(x»y)d(X|'yi> 

oH \\fi est une fonction de la forme 
' fviCx.y) = ai + Pix + Yiy (x, y) e e 

A Zvi(x.y) = i 

,Vi(x.y) = 0 (x,y)«e 

12. Precede selon I'une des revendications 1 a 11, 
caracterise en ce que la methode dif f erentielle 
d' optimisation des vecteurs deplacement est la methode 
de Gauss-Newton. 
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13. Precede selon I'une des revendications 1 & 11, 
caract^risS en ce que la mSthode dif f erentielle 
d« optimisation des vecteurs d^placement est 1' extension 

5 de Marquardt de la mSthode de Gauss-Newton. 

14. Proc^dS selon I'une des revendications pr6c6dentes, 
caract^ris6 en ce qu'une contrainte de compacitS est 
imposSe sur chaque element fini au moment de 

10 !• optimisation des vecteurs deplacement "du modele 
d' Pigments finis, ladite contrainte consistant a 
empecher la compacit€ de chaque element fini de tendre 
vers z6ro. 

15 15. Proc^dS selon la revendication 14, caract^ris6 en 
ce que la contrainte de compacitS sur un element fini e 
de sommets Pi, Pj, Pjc et de compacite C<Pi,Pj,Pk) est 
dSfinie par la relation suivante: 

C(Pi + dp. ,Pj + dp. ,Pk + dp^^ ) > K * C(Pi,Pj.Pk) 

20 oil dp..dp. etdp^ representent les vecteurs deplacement 
des sommits Pi, P j , Pk au cours de l'€tape 
d ' optimisation, 

et K est un parametre de compacitS. 
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16. Precede selon la revendication 14 ou 15, 
caracterisfe en ce que 1' optimisation des vecteurs 
deplacement sous contraintes sur les 61€ments finis est 
rgsolue par la technique dite du Lagrangien augmente. 

17. Precede selon la revendication 16, caract§ris4 en 
ce que les contraintes sont employees sous une forme 



linearisee dans la technique dite du Lagrangien 
augmente. 

18. Precede selon I'une des revendications 12, 13 ou 
17, caracterise en ce que les methodes d • optimisation 
des vecteurs deplacement utilisent une technique LDL*^ 
profil. 
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